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ABSTRACT 

Context. Gravitationally lensed quasars constitute an independent tool to derive the Hubble constant through the time-delays of their 
different images; they offer as well the opportunity to study the mass distribution and interstellar medium of their lensing galaxies 
and, through microlensing due to stars in their lensing galaxy, they also allow one to study details of the emitting source. 
Aims. For such studies, one needs to have an excellent knowledge of the close environment of the lensed images in order to model the 
lensing potential: this means, ideally, observational data over a large field-of-view and spectroscopy at high spatial resolution. 
Methods. We present VIMOS integral field observations of 27"x27" fields of view around four lensed quasars: HE 0230-2130, 
RX J091 1.4+0551, H 1413+117 and B 1359+154. Using the low, medium and high resolution modes, we study the quasar images 
and the quasar environments. Because of the complexity of the reduction of integral field datasets, we provide a detailed report of the 
data reduction. 

Results. Comparison between the quasar spectra of the different images reveals differences for HE 0230-2130, RX J09 11.4+0551 
and H 1413+117: flux ratios between the images of the same quasar are different when measured in the emission lines and in the 
continuum. We have also measured the redshifts of galaxies in the neighborhood of HE 0230-2130 and RX J09 11.4+0551 which 
possibly contribute to the total lensing potential, given their close proximity to the line-of-sight toward the quasars. 
Conclusions. A careful analysis reveals that microlensing is the most natural explanation for the (de)magnification of the continuum 
emitting region of the background sources. In HE 0230-2130, image D is likely to be affected by microlensing magnification; in 
RX 1091 1.4+0551, images Al and A3 are likely to be modified by microlensing de-magnification and inH 1413+117 at least image 
D is affected by microlensing. We have not been able to recover any microlensing information or galaxies close to the line-of-sight in 
the neighborhood of B 1359+154. 

Key words, gravitational lensing - techniques: spectroscopic - quasars: individual: HE 0230-2130, RX J091 1.4+0551, H 1413+1 17, 
B 1359+154 - cosmology: observations 



$_] 1 . Introduction Its study may allow us to statistically understand the flux ratios 

in order to use them to constrain the "macro" model of the lens. 

Gravitationally lensed quasars are a very useful tool to study Microlensing also provides valuable information about the 
the intergalactic medium, the mass distribution and substruc- background source object and the lens itself. By coupling pho- 
tures of the lensing g alaxies jgadeau et iFalcoetalJ tometnc or spectroscopic observations to theoretical rmcrolens 

1999; Wucknitzetal. 2003; YoneharlZS T5I1). They are ln 8 simulations (analytical or numencallone is able tode 

rive the size (e.g.. [W ambsgans s et al.l 11990; Yonehara 2001 



used in the quest for cosmological parameter s (Brada c et al 



2004; Meneghetti et al. 2005; Taylor 2005; Mortsell & Sunesson 



Kochanek et al. 2006) and shape (e.g., [Mineshige & Yonehara 
1999) of th e background source, a s well as the transverse ve- 



2006), in particular the Hubble constant (Ho) through time- 

delay measurements (lRefsdallll964l) . A list of m easured time- loclt y ( e -g- Gil-Merino et al. 2005J) and mean mass of the mi 

delays currently available has been established dOgurill2007i) . crolenses (e.g., |Kochanek|2004D . 

The Hubble constant can be derived using the time-delays be- The programme we are reporting upon in this paper was con- 

tween the images of the quasar and a model of the mass dis- ducted at ESO in period 74 (proposal IDs: 074.A-0152(A) and 

tribution of the lens, for which the amplification is directly 074.A-0152(B), PI: C. Faure) and intended to spectroscopically 

comparable to the observed flux ratio between the images, study in detail four gravitationally lensed quasars (HE 0230- 

However, lensing galaxies do not provide "smooth" potentials, 2130, RX J091 1.4+0551, H 1413+1 17 and B 1359+154), search 

but rather consist of stars. Therefore, there is often a supple- for the occurrence of microlensing and study their close neigh- 

mentary amplification caused by stars in the lensing galaxies borhood. 

close to the line-of-sight. This phenomenon is known as mi- Using the integral field spectroscopy mode of the Visible 

crolensing dChang & Refsdallll979t [W ambsga nss & Paczvnskil Multi Object Spectrograph (VIMOS) instrument installed on the 

1 199 It ISchechter & Wambsganss! l2002t ISchneider et al.l 12006). ESO/VLT-UT2 telescope in Cerro Paranal (Chile), we have ob- 
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tained 3D spectroscopy of the fields around the four targets. 
These are briefly introduced in Sect. [2] The dataset is presented 
in Sect. [3] The data reduction of the integral field spectroscopy, 
in its low, medium and high resolution modes, is detailed in Sect. 
|4] The identification and extraction of the object spectra in the 
quasar fields, together with their analysis, are reported upon in 
Sect. [5] Finally, the microlensing analysis is presented and dis- 
cussed in Sect. [6j followed by a summary and conclusions in 
Sect. |7] 

2. The targets 

The lensed quasar targets were selected based on the fact that 
the presence of a galaxy cluster/group was either known or 
suspected in the line-of-sight of the quasars. Furthermore, for 
RX J09 11. 4+0551 and H 1413+117, microlensin g activity was 
suspected dBade et al.lll997HAngonin et al.lll990l) . 

2.1. The quadruple quasar HE 0230-2130 

The quadruple lense d quasar HE 0230-2130 was discovered by 
Wisotzki et al. ( 1999) in the course of the Hamburg/ESO survey 
(Wisotzki et a 1. 1996). The quasar is at redshift z=2. 16. The sep- 
aration between the four quasar images ranges from Aab=0.74" 
to Aad=2.1". A galaxy ove rdensity was discovered 40" south- 
west of the quasar images dFaure et ai1l2004l) . lEigenbrod et al.l 
(2006) identified the redshift of the two main lensing galaxies 
that lie between the quasar images respectively at z^gi =0.523 + 
0.001 and z iG i =0.526 + 0.002. The images of the lensed quasar 
are displayed in the left insert of the right panel in Fig.Q] 

2.2. The quadruple quasar RX J091 1.4+0551 

The quasar RX J091 1.4+0551 w as originally selec ted from the 
ROSAT All-Sky Survey (R ASS) dBade et alJll995h . It was clas- 
sified bv lBade et al] d 19971) as a multiply imaged quasar with at 
least three images at redshift z=2.80. In an analysis with higher 
resolution imaging f rom the Nordic Opt ical Telescope (NOT, La 
Palma) observations. iBurud et al.l dl998l) discovered a fourth im- 
age of the quasar. The lensed quasar is described as three very 
close images (Al, A2 and A3, with A A i A2 =0.48", A A1A3 =0.96") 
and an isolated fourth image (B) located A / iib=3.05" west of 
image Al. The lensing system comprises an elongated lens- 
ing galaxy located at Aaig=0-86" and a galaxy cluster cen- 
tered at a distance of 38" from the quasar images. The redshift 
of the lensing galaxy was mea s ured from Keck ob servations: 
z/=0.769 (see iKneib etaT]|2000l) . Uneib et alj d2000t) character- 
ized the galaxy cluster located south-west of the quasar images at 
z=0. 769+0. 002, based on the spectroscopic measurements of 24 
galaxy members. The measured velocity dispersion of the galaxy 
cluster is of: <x= 836^00 km s _1 . The quasar images and the 
lensing galaxies are displayed in Fig. [2] 

2.3. The quadruple quasar H 1413+117 

The quadruply imaged quasar H 141 3+117 was first iden tified 
as a gravitationally lensed quasar by Mag ain et al.1 (Q988). It is 
one of the rare multiply imaged Broad Absorption Line (BAL) 
quasars. The quasar is at a redshift z=2.55, while the spectro- 
scopic redshift of the lensing galaxy is still unknown. The pres- 
ence of a gala xy cluster in the d irection of the quas ar was first 
suspec ted by IKneib etail dl998l) and confirmed by iFaure et al] 
(2004), at a photometric redshift z~ 0.8 + 0.3 (using up to six 



photometric bands). Additional absorption systems were identi- 
fied in the BAL qua sar spectrum at redshifts z=1.66 and z=1.44 
dMagain e t al. 1988). The different components of the system are 
displayed in Fig. [3] 

2.4. The multiple quasar B 1359+154 

The lensed quasar B 1359+15 4 was discovered a s part of the 
CLASS survey and discussed bv lMvers et al.l d 1999b . The quasar 
is at redshift z=3.24 an d has a ma x imum image separation of 
1.7". As identified by iRusin et al.l d200ll) . this system shows 
at least six images of a radio source and its star forming host 
galaxy. The main share of the lensing potential is produced by 
three lensing galaxies inside the Einstein ring with currently 
no spectroscopic redshift available; however, the three galaxies 
have been identified as the core of a galaxy group possibly at 
redshift z^ 1 based on optical colors and mas s estimates derived 
from models of the system (Rusin et al. 200lh. 

3. The VIMOS dataset 

This new set of observations was collected in service mode be- 
tween October 2004 and March 2005. The data were taken with 
the Integral Field Unit (IFU) using the low, medium and high 
spectral resolution modes. 

3.1. The VIMOS integral field unit 

The VIMOS integral field unit (IFU) consists of 6400 fibers cou- 
pled to microlenses, divided into four quadrants. Each quadrant 
is made of four sets of 400 fibers, called "pseudo-slits". Each 
pseudo-slit corresponds to a 20x20 pixels area within each quad- 
rant of the IFU head (80x80 pixels in total). 

The instrument can work in three resolution modes: 
low resolution (LR, spectral resolution=2 10-260, dis- 
persion=5.3 A pix~'), medium resolution (MR, spectral 
resolution=580-720, dispersion=2.5 A pix -1 ) and high 
resolution (HR, spectral resolution=2500-3100, disper- 
sion^. 6 A pix ). In the LR mode all fibers are used whereas 
in the MR and HR modes only one pseudo-slit is used per 
quadrant. Therefore, the size of the field-of-view depends on 
the spatial sampling. We have selected the sampling so as to 
cover a 27" x 27" field size (0.33"/fiber in LR and 0.67"/fiber in 
MR/HR). The targets were observed through the LR-blue, MR 
and/or HR-red/orange grismfl 

We have coupled the LR-blue grism with the OS -blue filter 
to cover a wavelength range of /l=3700-6700 A. The MR- and 
HR-red grism observations were coupled with the GG475 filter 
to cover a wavelength range of ,1=5000- 10000 A and ,1=6350- 
8600 A respectively, while the HR-orange grism was coupled 
with the GG435 filter to cover a wavelength range of ,1=5250- 
7550 A. 

3.2. The observing runs 

A summary of the observing runs used for this analysis is dis- 
played in Table 1. For RX J091 1.4+0551 and HE 0230-2130 
we have retained only the exposures conducted under a seeing 

1 For HR, the observations through quadrants 1, 2 and 3 were ob- 
tained with the red filter, while there was no red grism for quadrant 4. 
Here we used the orange grism instead, see VIMOS Handbook 30 June 
2005. 
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Table 1. Observation log (all observations performed between 
October 2004 and March 2005). 



Object 


Mode 


Night 


FWHM 

" 


Ait- 
Mass 


Exp 

sec 


HE 0230-2130 


LR Blue 


Nov. 17 


0.60 


1.24 


2355 






Nov. 17 


0.58 


1.50 


2355 






Dec. 5 


1.08 


1.02 


2355 


HE 0230-2130 


HR Red 


Oct. 16 


0.60 


1.18 


2205 






Oct. 16 


0.72 


1.07 


2205 






Oct. 16 


0.87 


1.02 


2205 






Oct. 16 


0.61 


1.09 


2205 






Nov. 13 


0.54 


1.24 


2205 


RX J0911.4+0551 


LR Blue 


Dec. 8 


0.58 


1.28 


2355 






Dec. 9 


0.67 


1.31 


2355 






Dec. 15 


0.83 


1.29 


2355 


RX J0911.4+0551 


HR Red 


Dec. 5 


0.47 


1.25 


2205 






Dec. 5 


0.46 


1.41 


2205 






Dec. 10 


0.52 


1.24 


2205 






Dec. 10 


0.50 


1.40 


2205 


H 1413+117 


LR Blue 


Mar. 18 


1.27 


1.25 


2355 






Mar. 18 


1.13 


1.26 


2355 






Mar. 18 


1.21 


1.64 


2355 


H 1413+117 


MR Orange 


Mar. 17 


0.83 


1.24 


2205 






Mar. 17 


1.22 


1.29 


2205 






Mar. 17 


1.03 


1.42 


2205 






Mar. 17 


1.12 


1.68 


2205 


B 1359+154 


MR Orange 


Feb .13 


1.40 


1.24 


2205 






Mar. 4 


0.81 


1.33 


2205 






Mar. 15 


0.59 


1.30 


2205 






Mar. 15 


0.39 


1.33 


2205 






Mar. 16 


0.88 


1.38 


2205 



below 1.1" and an airmass less than 1.7. For H 1413+117 and 
B 1359+154 we have considered the observations including at 
least three quadrants and obtained with an airmass below 1 .7, as 
well. This represents total exposure times of 7 ks in LR and 1 1 
ks in HR for HE 0230-2130, of 7 ks in LR and 9 ks in HR for 
RX J09 1 1 .4+055 1 , of 7 ks in LR and 9 ks in MR for H 1 4 1 3 + 1 1 7 
and of 1 Iks in MR for B 1359+154. 



4. The data reduction 

The data reduction has been performed mainly using the 
VIMOS Interactive Pipeline & Graphical Interface (VIPGI 
IScodeggio et alJl2005l:IZanichelli et al.ll2005l) . While most infor- 
mation concerning the data processin g can be retrieved from the 
VIPGI Cookbook (Garilli et al. 2007) and some important issues 
are discussed in ICovone et alj d2006l) . we find it necessary to 
provide additional details on some of the most complex reduc- 
tion steps in order to make the VIMOS IFU data reduction more 
transparent for future VIPGI users. 

4.1. First adjustments 

The location of the object spectra and the wavelength calibra- 
tion both require an accurate description of the instrumental dis- 
tortions. The pattern of these distortions is provided in the data 
header. However, the instrument flexures change with the point- 
ing direction and therefore these values can only be used as 
"first guesses" for a further improvement by polynomial fits. The 
VIPGI software offers the opportunity to tailor these first guesses 
separately for each frame in the observation run. 

During this step, individual fibers can be flagged as "dead" 
if necessary (i.e. zero fiber transmission/efficiency). The flagged 
fibers will not be considered in any of the subsequent fittings. 
If the flagging is not made properly, it yields a distorted image 
of the astronomical object as illustrated for example in Fig. [4] A 
quality control routine allows us to check the adjustments in the 



Fig. 4. Resulting images of a point-like object in the HR mode 
for incorrect (left panel) and correct (right panel) flagging of the 
fibers. In the left panel dead fibers have not been flagged, hence 
during the first guess adjustment, the software assigns to these 
fibers the flux of the adjacent fiber. 



different science exposures. The corrected first guesses are then 
applied to the entire scientific and calibration dataset. 

In the second quadrant of the LR mode, the first 80 fibers 
were flagged as dead because of their low transmission in all 
frames. 

4.2. Spectra tracing and extraction 

After the first guess correction, one needs to trace the spectra 
in the different exposures. To do so, we use a high signal flat 
field exposure obtained immediately after the science exposure. 
In the LR mode, we use the global fit parameters proposed in 
the image headers, as they allow a good tracing; for the pixel 
binning we use 20 pixels to trace the spectra from their start 
to their end. For the MR/HR frames, better results are obtained 
ignoring the global fit distortion model: the IFU binning then 
needs to be larger than 20 pixels (the best fit was obtained for 
~50 pixels). From this fit, a table of spectra positions is built. 
The position table has to be visually inspected to check whether 
the spectra are correctly traced. 

The inverse dispersion solution is calculated next in the "cre- 
ate master lamp" routine. This is done by fitting a third order 
polynomial to the lamp frame, starting from the first guesses. 
The inverse dispersion solution was checked visually over the 
raw frames as well as through the distribution of the root mean 
square (rms) residuals of the fit. The rms residuals are of about 
0.7 A ± 0.7 A for the LR frames and 0.09 A ± 0.06 A for the 
MR/HR frames. 

Independently, the science frames are corrected by a master- 
bias and over-scan trimming. Then, the object spectra are ex- 
tracted from the science exposures. At this stage, cosmic rays 
are detected and eliminated and the wavelength calibration is 
applied. 

4.3. Last calibration and exposure combination 

The final phase of the reduction process includes many steps. 
First, the software applies a fiber-to-fiber correction in which the 
transmission of each fiber is calibrated. This fiber-to-fiber cor- 
rection is calculated using a bright and isolated sky-line present 
in all fibers. The transmission is corrected so that each fiber has 
the same flux all along the sky-line. To perform this correction 
we have selected the 5577 A sky-line in LR-blue, MR and HR- 
orange and the 6300 A skyline in HR-red. 

The sky contributi on is dealt with in a statistical way 
dScodeggio et alj r2005l) by subtracting the mode spectrum in 
each pseudo-slit. This gives satisfactory results if at least 50% 
of the fibers in the frames relate to pure sky, which is the case 
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for the current dataset. The spectro-photometric calibration of 
the science frames was made using the observation of standard 
stars, which accurately corrects for the relative flux. 

The information available in the science frame headers con- 
cerning the position of the pointing for our data is not accurate 
enough to be used as a guide for combining the dithered expo- 
sures. Therefore, we have computed the shifts between the dif- 
ferent exposures using the brightest objects in the frame. This is 
a delicate step because of a Point Spread Function (PSF) shape 
problem which is discussed in Sect. 4.4. 

Then, the different exposures are combined using a median 
filter to eliminate the dead pixels and remaining cosmic rays. 
The resulting product is a datacube: for each pixel within the 
27"x27" region of the sky, a spectrum is available. During this 
combination, a fringing correction may be applied. 

The datacube is then corrected for differential atmospheric 
refraction. Applying a theoretical differential atmospheric re- 
fraction correction was attempted but does not give satisfactory 
results. Instead, we have computed and applied an empirical dif- 
ferential atmospheric correction: in most cases, the quasars are 
very bright in the datacube, so we have been able to fit a poly- 
nomial to the highest flux position of the quasar image for each 
wavelength. We then knew the position correction to be applied 
to each wavelength cut of the datacube. 

Images of the fields around HE 0230-2130, 
RX J091 1.4+0551 and H 1413 + 117 are shown in Figs [TJ 
[2] and [3] The images were created by adding all the monochro- 
matic slices in the datacube, excluding those corresponding to 
strong sky emission lines, and excluding the edges of the grisms. 
In the case of B 1359+154 we only have MR observations with 
one missing quadrant and low signal-to-noise ratio (S/N ~4), 
hence we could not create such an image. 

4.4. An additional PSF problem with VIMOS 

After the reduction steps and the corresponding checks described 
above, the PSF in the science target and standard star datacubes 
still appeared elongated along the East-West direction for all in- 
dividual exposures. We could not find any VIMOS documenta- 
tion reporting on this problem; yet, it has been confirmed by ESO 
staff (private communication) that as the IFU head is placed at 
the edge of the VIMOS field-of-view, the image of a point source 
is subject to an aberration and appears slightly elongated in one 
direction (see for example image B of RX J091 1.4+0551 in the 
left panel of Fig. 2). This seems to be an identified problem to 
be mentioned in the VIMOS documentation. Besides the time 
investment in understanding the reason for the elongation, it is 
very difficult to quantify how much this issue impairs the im- 
age resolution in the datacubes, as this is a combination of other 
effects as well (seeing, signal-to-noise and dithering). 

However, and as shown below, it is still possible to extract 
relevant spatial information from the data. 

5. Identification of the quasar and galaxy spectra 

A number of tools were used to visualize and analyze VIMOS 
datacubes. We used mainly the Euro3D software (Sanchez 2005, 
http://www.aip.de/Euro3D) that allows us to inspect the dat- 
acube, both in an interactive and an automated manner (guidance 
on the use of the software can be found in the Euro3D documen- 
tation). 

We have used existing deep images of the target field- 
of-views to identify and locate objects that are faint in the 



VIMOS rebuilt images (27"x27"). We have used VLT/FORS1 
R band data to locate objects in the field around HE 0230- 
2130, HST/WFPC-F555w and I814w for RX J091 1.4+0551 and 
HST/WFPC-I814w, ISAAC J and ISAAC K for H 1413+117 
(HST programmes Pis, respectively: J. Hjorth and A. Westphal 
and VET programme PI: C. Faure). 

The objects for which we have detected a signal are 
shown and labeled in Figs. QJ and for HE 0230-2130, 
RX J091 1.4+0551 and H 1413+117 respectively. For redshift 
analysis we use the e^I software coupled with spectral templates 
dKinnev et al.| [l996). plus iraf routines. We provide below a de- 
tailed analysis of the different spectra. 

5.1. The field of HE 0230-2130 

5.1 .1 . The quasar spectra 

The spectra of the quasar images have a very high signal-to- 
noise ratio (S/N ~40 in the continuum). However, because of 
PSF shape problems and of the small angular separation between 
the quasar images, it is not possible to separate all the spectra 
of the different quasar images. We have partitioned the images 
in two groups, A+B (from now on AB) and C+D (from now 
on CD), which are visually resolved in the final datacube (see 
Fig.QJfor labels). Based on the Lyman or(1215A) + NV(1240A), 
SiIV+OIV](1398A), CIV(1550A) and CIII](1908A) emission 
lines, the spectra from both AB and CD provide a source red - 
shift z=2.163 + 0.003, in agreement wifh lWisotzki et all dl999h . 
In the two spectra an absorption system is detected at the same 
redshift as the quasar, likely originating from the quasar host 
galaxy. No other obvious absorption system is detected in the 
spectra (see Fig. [5]). 

5.1.2. Galaxies in the field around HE 0230-2130 

We have extracted spectra at the location of the main lensing 
galaxies (LG1 and LG2 as shown in Fig. Q} that are highly 
contaminated by quasar emission. Therefore, a scaled version 
of the spectrum of the isolated quasar image C has been sub- 
tracted. Doing so, we are able to detect absorption lines from 
both lensing galaxies. We measure the redshift of LG1 at 
z=0.521 +0.004 and of LG2 at z=0.524 ±0.003. These redshifts 
are in agreement with those obtained by Eigenbr od et alj (2006) 
(zlgi =0.523+0.001 ; zlg2=0.526±0.002). 

Another galaxy is seen in the field, north west of the quasar 
image A (Aa=-15", A<5=9.2", Galaxy.l in Fig.QJi. Its spectrum 
shows the Call K, Ca II H and G-band absorptions, and also the 
[OIII] and [Oil] emission lines (see Fig. [6j. From these lines 
we measure a redshift of z=0.518 + 0.002, thus, this galaxy is 
probably a member of the group or cluster to which the lensing 
galaxies LG1 and LG2 belong. 

5.2. The field of RX J091 1 .4+0551 
5.2.1 . The quasar spectra 

Because of the small separations between images Al, A2 and A3 
(A A ia2 =0.48", A A ia3 =0.96", A A 2A3 =0.61") and the distorted 
shape of the PSF, it was not possible to separate the spectra of 
the three quasar images. Therefore they are treated as a single 
object called A (see Fig. [2] left panel). 



2 Redshift analysis software by the authors of VIPGI: 
http://cosmos.iasf-milano.inaf.it/pandora/EZ.html 



VIM OS 



% 

\ 



CD 



AB 



T. Anguita et al.: Integral field spectroscopy of lensed quasars 

FORS1 



HE 0230-2130 



c ■ 



LCI 

b 



Galaxy_l 

® 



HE 0230-2130 



Fig. 1. The 27"x27" field-of-view around the multiple quasar HE 0230-2130: VIMOS IFU LR data cube (left panel) and FORS1 R 
band image (right panel). The components of the system are labeled as in the text. Orientation and scale are identical in the left and 
right panels. A FORS1 R band zoom of the system is displayed in the bottom left hand corner of the right panel. 
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Fig. 2. The 27"x27" field-of-view around the multiple quasar RX J091 1.4+0551: VIMOS IFU LR data (left panel) and WFPC 
F555w band image (right panel). The components of the system are labeled as in the text. Orientation and scale are identical in the 
left and right panels. A WFPC F814w band zoom of the system is displayed in the bottom left hand corner of the right panel. 




Fig. 3. The 27"x27" field-of-view around the multiple quasar H 1413+117: VIMOS IFU LR data (left panel) and ISAAC J band 
image (right panel). The components of the system are labeled as in the text. Orientation and scale are identical in the left and right 
panels. A WFPC F814w band zoom of the system is displayed in the bottom left hand corner of the right panel. 



In the LR-blue mode, the A and B quasar spectra exhibit shift z=2. 790 ± 0.009. As already mentioned by iBade et a"D 
Lyman c+NV+OV], SilV+OIV] and CIV emission lines at red- (1 19971) . we can also detect several absorption lines in both 
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Fig. 5. Low resolution spectra for HE 0230-2130. Top panel: 
combined spectra for images AB. Middle Panel: combined spec- 
tra for images CD. Bottom panel: flux ratio between the two 
spectra. Shaded regions indicate the spectral ranges dominated 
by intense sky line residuals. 



HE 0230-2130: Galaxy_1 (z = 0.518) 




3500 4000 4500 5000 5500 6000 6500 
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Fig. 6. Neighborhood of HE 0230-2130: Low resolution spec- 
trum of the galaxy in the field around HE 0230-2130 (labeled 
Galaxy.l in Fig. 2) and a scaled version of a template elliptical 
galaxy spectrum. From the lines shown we measure a redshift of 
z=0.518 + 0.002. 



Fig. 7. Low resolution spectra for RX J09 1 1 .4+055 1 . Top panel: 
summed spectrum for image A. Middle panel: spectrum for im- 
age B. Bottom panel: flux ratio between the spectra A and B. 
Shaded regions indicate spectral ranges dominated by sky resid- 
uals. 

Table 2. Values obtained for the five parameter Gaussian fits of 
the absorption lines detected in the spectra of RX J09 1 1 .4+055 1 . 
An asterisk (*) denotes the lines close to the 5577 A sky line. 



Redshift 


Line 


Center [A] 


FWHM [A] 


cr kms 1 


z=2.80 


Lya 










Civ 


5840.5 


15 


330 




Cm 


7185.5 


19 


340 


z=2.63 


Lya 


4404 


30 


940 




Civ* 


5618.5 


19 


430 




Cm 








z=2.54 


Lya 










C,v* 


5487.5 


28 


640 




Cm 








z=2.45 


Lya 


4163 


29 


900 




Civ 










Cm 









the A and B spectra; CIV and Lyman a at z=2.42 (or ~ 
-29,000 km s 1 with respect to the quasar), CIV, SilV, and 
Lyman a at z=2.63 (or ~ -13,000 km s _1 with respect to the 
quasar), and a few more absorption lines at wavelengths be- 
low 4500A. An absorption line at 5488A, if identified as CIV, 
would imply a redshift of z=2.54 (or ~ -22,000 km s _1 with 
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Fig. 8. High resolution red spectra for RX J091 1.4+0551. Top 
panel: total spectrum for image A. Middle panel: spectrum for 
image B. Bottom panel: flux ratio between the two spectra. Sky 
line residuals dominate the region redder than 7500A. 



RX J091 1 .4+0551 : Galaxy_1 (z = 0.484) 
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Fig. 9. Neighborhood of RX J091 1.4+0551: Spectrum of 
Galaxy_l and of an Sa spectral template for comparison. 
Identified lines are highlighted by dotted straight lines. They 
yield a redshift z=0.484 ± 0.001 



5.2.2. Galaxies in the field around RX J091 1 .4+0551 

We are unable to detect any absorption features from the lensing 
galaxy because of the intense contribution of the three quasar 
images. 

We identify three objects in the field: Galaxy _1, Galaxy _2 
and Spiral as labeled in Fig. [2] We are only able to measure 
the redshift of Galaxy_l for which we obtain a value z=0.484 
+ 0.001 (see Fi g. |9), indicating th at it is not a member of the 
cluster found bv lKneib et all d2000h at z=0.769. 

In the HR-red mode, we identify the galaxy labeled 
KCH2000-3 in Fig . [2] which had been previously detected by 
iKneib et aU yOOO) and measured at z=0.7618. Our determi- 
nation of the redshift gives z=0.762, in good agreement. This 
galaxy is part of the galaxy group discovered south-west of the 
quasar. 



respect to the quasar). It might correspond to the line identified 
bv lBade etail d 19971) at 5546A (z=2.57), although its proximity 
to the 5577 A intense sky line makes its redshift determination 
uncertain. The absorption lines have a velocity dispersions of 
<x < 1000 km s _I (see Table d, hence, the quasa r was classified 
as a mini-BAL quasar (e.g., Chartas et al. 2001). The physical 
origin of the absorption lines is still under discussion; the most 
likely interpretation is that these lines are due t o high-velocity 
outflo ws related to the quasar central engine (e.g. jReichard et all 
120031) . 

In the HR-red spectra of the A and B quasar images, we de- 
tect the CIII] emissio n line at the quas ar redshift, not seen in 
the previous spectra of lBade et all d!997l) . Thanks to the VIMOS 
high quality data, we also detect a few absorption features on top 
of the CIV and CIII] emission lines. Their FWHM are smaller 
than those of the other absorption lines previously indicated (see 
Table [2]). The narrow features could either be produced by out- 
flowing material from the quasar central engine, or signatures of 
the quasar host galaxy. 

In Figs. [7] and [8] the spectra of images A and B are shown 
in low and high resolution modes respectively, as well as the 
flux ratios between the spectra. Residuals in the flux ratios are 
discussed in Sect. 6. 



5.3. The field around H 1413+117 

The quasar images cannot be visually resolved in the two di- 
mensional cuts, either in the LR or in the MR spectroscopy. The 
low and medium resolution summed spectra for all four quasar 
images show the Lyman a, NV, SilV emission lines coupled 
with their broad absorption counterparts and the CIII] line. From 
the emission lines we measure the quasar redshift at z=2.554 ± 
0.002. 

The higher spectral resolution and longer wavelength range 
available in the medium resolution mode allows us to detect th e 
absorption systems already described by Mag ain et al.l (Q988). 
Taking advantage of the improved resolution and small width 
of the absorption lines we have been able to obtain very tight 
values for their redshifts. In the first absorption system, the Fell 
triplet (-5750A) and the Mgll doublet (-6850A) yield a red- 
shift z= 1.4381 ±0.0002; in the second system, the A1III doublet 
(-4950A), Fell triplet (-6300A), the Fell doublet (-6900A) and 
the Mgll doublet (-7450A) lead to a redshift z=1.657±0.001. 
Two other doublets are detected at ~6450A and ~6600A; they 
do not seem to be related and it is hazardous to assign them a 
specific identification. Low and mid-resolution spectra are dis- 
played in Fig. [10] 



8 



T. Anguita et al.: Integral field spectroscopy of lensed quasars 




Rest Vuivekmot.h [Al 

iaoo 2000 



H1413+117- 
z = 2.55 




5500 6000 6500 7000 

Wavelength [A] 



7500 8000 



Fig. 10. Top panel: LR-blue spectrum for H 1413+117. The 
spectrum was reconstructed by adding the contributions from 
the four quasar images. The dashed areas highlight regions af- 
fected by sky subtraction residuals. Bottom panel: MR-orange 
spectrum for H 1413+117 and the summed contributions from 
the four images. The lines identified with the z=1.44 and z=1.66 
absorption systems are indicated. Two other doublet systems are 
seen at -6450A and -6600A 



5.4. The field of B 1359+154 

For the quasar system B 1359+154 only the MR-mode was 
used. It being a faint and high redshift quasar, the signal-to-noise 
ratio is lower than in the other datasets (S/N ~ 4 in the contin- 
uum). 

However, we can recover a spectrum of the quasar by 
summing over the spatial elements in the two dimensional 
wavelength cuts (spaxels) of the different quasar images. The 
spectrum (displayed in Fig. PTT1) exhibits the Lyman a, NV, 
SilV+OIV], CIV and Hell lines, yielding a redshift z=3.235 + 
0.003. 

No absorption system is seen in the spectrum. In addition 
no galaxy is detected in the field. As no spectroscopic signature 
could be recovered from individual images, no further analysis 
was possible for this quasar. 



6. Interpretation of the spectral differences between 
the quasar images 

6.1. Spectral differences in HE 0230-2130 

In order to identify and analyze any differences between the 
spectra from images AB and CD, we first simply measure the 




Fig. 11. MR-orange spectrum of B 1359+154 (all spaxels with 
signal have been summed). Shaded areas correspond to regions 
affected by sky subtraction residuals. 



ratio between the two spectra. The ratio varies with wavelength, 
in particular it is different if measured over the continuum or 
over the emission lines. Furthermore, a mild slope in the con- 
tinuum flux ratio is observed. Notice that the edge of the grism 
(<4000A) is an untrustable region and has been dropped in the 
analysis. 

To study this phenomenon further, we developed a specific 
procedure to extract the spectra and spectral energy distribution 
(SED) of the four individual quasar images: a surface with four 
Gaussian profiles was fitted to each 2D monochromatic slice of 
the LR-blue datacube. Distances between the Gaussian centroids 
correspond to the distances between the quasar images as pro- 
vided in the CASTLES web pagfl from HST data. A further 
constraint is that the FWHM (in both directions: x and y) should 
be identical for the four profiles at each wavelength cut. The fit 
was performed using Levenberg-Markwardt routines in an itera- 
tive fashion (MPFIT, Craig MarkwardQ). In Fig.[T2]we show the 
extracted spectra for images A and D of HE 0230-2130, while 
their ratio is displayed in Fig. 13 (also reported in Fig. 14). The 
spectra of the isolated images A and D display emission line dif- 
ferences as well as a mild continuum slope difference. This is 
not surprising as image D of the quasar is located almost behind 
one of the lensing galaxies: LG1 (Fig. 2), therefore it is prone 
to be affected by individual stars and interstellar matter in LG1. 
Indeed, the flux ratios between images A and B and between im- 
ages A and C (not shown here) show little or no variation above 
the noise, contrary to image D. 

Intrinsic quasar variability is a possible, but unlikely, expla- 
nation for this phenomenon. Indeed, images A and D are sepa- 
rated by a time-delay and therefore any intrinsic variation of the 
quasar flux would be reflected at different epochs in images A 
and D. However, intrinsic variability of the quasar itself would 
more likely show selective emission line differences rather than 
similar magnification in all the emission lines. Indeed, differ- 
ent lines being emitted from differe nt regions in the q uasar are 
not likely to vary at the same time dKaspi et al.ll2000l) . Hence a 
more plausible explanation for the emission line differences is 
microlensing by stars in the lensing galaxy, affecting preferen- 
tially image D of the quasar. 

Microlensing is expected to affect the continuum emission 
more strongly than the broad line emission, because the con- 
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Fig. 12. Spectra of the individual images A and D of HE 0230- 
2130. The normalized spectrum for each image is displayed. The 
spectrum for image A has been shifted upward by 1 .5 flux units 
for comparison. Shaded areas highlight regions affected by sky 
subtraction residuals. 

tinuum emitting region is thought to have a size comparable to 
those of the Einstein rings of stellar microlenses, whereas the 
broad line region is much larger (e.g JSchneider & Wa mbsganss 
1990; AbajasetaD [2002). To demonstrate this: based on the 
virial theore m, the size of the br oad line region (/?blr) is esti- 
mated to be dPeterson et al.lll998l) : 



* BLR = °- 57pc llowr) 



Mu 



(. 

\ 1O 8 M 



(1) 



where vfwhm is the velocity derived from the FWHM of an 
emission line and M BH is the black hole mass. Assuming a rea- 
sonable size for the continuum emitting region (R cont ), of 1 , 000 
Schwarzschild radii (r s ), we get: 



/ M] 



* cont = 9.1X10 P c^— 



(2) 



It is obvious that R cont <k ^blr, independently of the black 
hole mass. In this case, the continuum of image D could be mag- 
nified by microlensing, making the emission lines "sink" in the 
continuum and, therefore, making them appear small in compar- 
ison to the emission lines of image A. 

Could microlensing also explain the slope difference (chro- 
matic effect) seen between the two images A and D? Accretion 
disk models predict that continuum emissi on from the central 
part i s bluer than that from the outer part (Shakura & Sunvaevl 
1 1973b . In general, the innermost part of an accretion disk 
(blue-emitting) is more strongly affected by microlensing than 
the outer part (red-emitting) as the Einstein ring of mi- 
crolenses matches better the size of the inner accre tion disk 
dWambsganss & Paczvnskil 1 1 99 it lYonehara et~ai1ll998l) : there- 
fore, to first order and in a naive vision, an image affected by 
microlensing should exhibit a bluer continuum. In image D, the 
situation is reversed: the quasar image thought to be affected 
by microlensing (image D) appears redder than the image not 
affected by microlensing (image A). Let us examine how this 
could happen. 

In Fig. 13 we illustrate how microlensing can change the 
slope of the flux ratio between images A and D, by using a sim- 
ple model for the caustic. We adopt a sta ndard accretion disk 
model for the continuum emitting region (Sha kura & Sunvaevl 
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Fig. 13. Top panel: Schematic of the straight line caustic with re- 
spect to the source, d is the distance measured between the center 
of the source and the caustic. Bottom figure: flux ratio between 
HE 0230-2130 images A and D. Superimposed on the flux ratio 
observed curve, the straight lines are the modeled flux ratios in 
different situations, that is for different distances d. The solid line 
is the modeled flux ratio without microlensing (i.e. d=oo). The 
dotted, the dashed, and the long-dashed lines correspond, respec- 
tively, to the modeled flux ratios with microlensing for distances 
between the caustic and the source center as described in the 
legend. The shaded area denotes the blue extremum of the grism 
(untrustable region). 



119731) . The black hole mass is set at 1O 8 M , the accretion rate is 
set fixed at the critical accretion rate and the inner and outer radii 
are set at 3r s and 1000r s respectively. Additionally, we assume a 
single straight line caustic and appl y the following approx imated 
formula for magnification (jj) (e.g., Schneid er et al.ll 1 992b : 



-1/2 



-X-scale J 

yu a (x<0) 



+ /j. d (x>0) 



(3) 



where jc sca i e is the scale length of the caustic and is chosen to be 
identical to the Einstein ring radius of a lensing object with mass 
1M Q . The parameter /u a represents the total magnification except 
for a pair of micro-images which appear/disappear at the caustic 
crossing. For image A and image D, the /i a values are arbitrarily 
set to be 3.0 and 2.0 in order to match the observed continuum 
flux ratio ~1.5 (see Fig. 13). 

Microlensing is not the only possible explanation for the 
chromatic effect. Galactic dust extinction through the lensing 
galaxy is also an alternative. In this case, galactic dust would 
block the blue light from image D, making this image look 
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Fig. 14. A normal extinction law, Ry = 3.1, for dust in the Milky 
Way (solid line) and an extinction law as in the Small Magellanic 
Cloud (dotted line). The shaded area denotes the blue extremum 
of the grism (untrustable region). 

redder in comparison to image A, and resulting in a decreas- 
ing flux ratio towards the redder part of the spectrum, as ob- 
served between images A and D. We evaluate the chromatic 
feature w hich is to be expec ted from differential dust extinc- 
tion (e.g.. lJean & Su rdei 1998). For these calculations, we con- 
sider a normal extinction law, Ry = 3.1, for dust in the Milky 
Way (Ca rdelli et al.lll989l) and an extinc tion law as in the Small 
Magellanic Cloud dGordon et al.ll2003|). In both c ases, we as- 
sume a local gas-to-dust ratio ( Bohlin et al.lll978h Ay = 5.3 X 
\0" 22 Nh, where Ay and Nh are, respectively, the dust extinction 
in the V band (in magnitudes) and the column density of hydro- 
gen gas in cirT 2 . We also assume that the difference in the gas 
column density in front of the images is 5 x 10 20 cirT 2 in the 
case of the Milky Way, and is 5 x 10 21 cirT 2 in the case of Small 
Magellanic Cloud. Fig. [14] shows the extinction expected from 
these two models. Both of them seem to trace at least part of the 
slope that we see in the continuum flux ratio and therefore galac- 
tic extinction is a plausible explanation for the phenomenon. 

Thus, in the case of HE 0230-2130, our preferred interpreta- 
tion for the difference between the spectra of images A and D is 
microlensing due to the lensing galaxy LG1 and affecting image 
D, as it explains in the spectrum of image D, both its emission 
line dimming and its continuum reddening. However, this is not 
the only interpretation: in particular the effect of the interstellar 
matter and dust in the lensing galaxy LG1 should not be over- 
looked. 



6.2. Spectral differences in RX J0911. 4+0551 

The different images of quasar RX J09 11.4+0551 also exhibit 
spectral flux differences. However due to the small separation of 
the images, it is impossible to apply the same extraction proce- 
dure as for HE 0230-2130. In order to quantify the differences 
between the spectra, we have performed Gaussian profile fits to 
the three main emission lines in the LR-blue wavelength range. 
Using the continuum levels obtained from the Gaussian fits we 
have derived the equivalent widths (EW) of the emission lines, 
and considered them as a measure of their strength. This was 
performed using: 



EW 



-f 



Table 3. Parameters of the Gaussian fits for the emission lines 
in the system RX J09 1 1 .4+055 1 : amplitude and continuum level 
are expressed in 10~ 15 erg position centers, Gaussian cr and 
equivalent widths in A. The equivalent widths are measured in 
the rest frame. 



LINE 




a 




Sin 




Civ 


IMAGE 


A 


B 


A 


B 


A 


B 


Amp. 


4.42 


0.59 


0.72 


0.10 


1.86 


0.24 


*Cmt [A] 


4639 


4636 


5292 


5287 


5853 


5843 


<r[A] 


75.95 


73.69 


33.03 


31.63 


55.61 


59.35 


Cont. 


2.85 


0.47 


2.94 


0.50 


2.70 


0.46 


EW [A] 


81.83 


63.69 


5.43 


4.50 


24.93 


20.13 



where S (A) is the flux at wavelength A and I is the continuum 
level for the chosen wavelength range (in this case the range of 
the emission lines, delimited by A\ and /I2). The parameters are 
displayed in Table 3. 

The differences between images A and B are confirmed by look- 
ing at the equivalent width values which are systematically larger 
in image A. While the continuum ratio is around 5.9, the emis- 
sion line ratio goes up to 7.7. Emission line residuals are seen 
in the flux ratio between images A and B, but in contrast to 
HE 0230-2130, the emission lines appear to be stronger in the 
image that is closer to the lensing galaxies. As argued in Sect. 
6.1, we think that microlensing is the most probable explanation 
for the emission line residuals, but in this case, given the dis- 
tances of the images to the lensing galaxy, we are rather dealing 
with microlensing demagnification in the image complex A. 

To probe this interpretation, we have built a "macro" lens 
model for RX J09 1 1 .4+055 1 , using an e lliptical effective lensing 
potential (Blandford & Kocha nekll 19871) with external shear. 

Using this lens model, we have evaluated the effective shear 
(jeff) and the convergence (/c* f ^) on all the images of the 
quasar, in order to characterize the properties of the Fermat 
travel-time surface dBlandford & Narayan 1986): minima, max- 
ima or saddle points. The results are shown in Fig. 15. Since im- 
ages Al and A3 are located at saddle points, they are more likely 
to be demagnified by microlensing dSchechter & Wambsganssl 
2002), therefore it is quite plausible that the continuum emitting 
region of image complex A (Al, A2, A3) is overall demagnified. 

Of course, the effective shear and convergence depend on the 
convergence provided by a smooth distribution of matter (k c ) and 
a clumpy stellar component (k*) in the following way: 



eff 
K' - 



Jeff = 



1 - K c 

J 
1 - K c 



(5) 



(S(A)-I) 



dA 



(4) 



However, the quasar images Al and A3 lie at saddle points, 
no matter which fraction of k c and a:* is incorporated in the lens 
model (see Fig. [TBI. 

Regarding the slope of the continuum flux ratio for 
RX J091 1.4+0551, a mild tilt is also observed. We have mea- 
sured the slope of this variation by adjusting a straight line over 
the whole spectral range (see Fig. [T6t. The fit has been adjusted 
over the continuum ratio which was obtained after subtracting 
the Gaussian fits of the emission lines from the total flux ratio. 
We have measured a slope of (-1.75 ± 0.15) x 10~ 4 A~'. Due to 
the evidence for microlensing demagnification, as argued above 
to explain the emission line residuals, we are led to consider that 
microlensing is also responsible for the chromatic distortions, 
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Fig. 15. Results of a simple lens model for RX J091 1.4+0551: 
Different values for the effective shear versus effective conver- 
gence. The symbols correspond to different images of the quasar: 
black squares: Al, black triangles: A2, crosses: A3 and stars: B. 
From left to right, the contribution of k c to varies from 100% 
and 0%. All combinations place the images Al and A3 in the 
zone of saddle points. 
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Fig. 17. Individual spectra for the four images of the quasar 
H 1413+117 after the fitting procedure. They have been normal- 
ized and shifted vertically for easier comparison (by 2.5 with re- 
spect to the bottom one). Shaded areas highlight regions affected 
by sky residuals. 
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Fig. 18. The normalized spectrum of image D of H 1413+117 
(thin line) and the mean spectrum of images A, B and C (thick 
line). Shaded areas show regions affected by sky line residuals 
may. The difference in the relative amplification between the 
continuum and the emission lines is evident, especially at the 
SilV and CIV lines. 



Fig. 16. A linear fit to the continuum ratio for the lensing system 
RX J09 1 1 .4+055 1 . A mild slope can be seen. 

making the complex image A appear bluer. However, this mea- 
sured effect is of small significance and impossible to disentan- 
gle from differential extinction with the present dataset. 

6.3. Spectral differences in H 1413+117 

In order to extract an individual spectrum for each image of the 
quasar we apply the procedure described earlier in Sect 6.1. As 
before, owing to the small separation of the images, it is possible 
to do so only with data obtained in the low resolution configura- 
tion. The spectra of the four images are displayed in Fig. [171 

We notice some differences between the spectra. The quasar 
image D exhibits the most obvious differences when compared 
to the other images (see also FigfTSI): the CIV line clearly shows 
a smaller intensity with respect to the continuum, while the 
SilV+OIV] line has practically vanished. This again is consis- 



tent with the continuum of image D being magnified by mi- 
crolensing. 

iKavser et al] dl990l) suggested microlensing as an explana- 
tion for photometric variations seen in image D of the quasar. 
Using integral field spectroscopy observations, lAngonin et al] 
( I 1990 s ) found spectral differences between image D and the other 
three images. They suggested that the differences were produced 
by broad absorption line clouds being selectively magnified by 
microlensing. A co mparison between the p resent VIMOS spec- 
tral dataset and the Angoni net al] ([1990) dataset, in particu- 
lar for image D, shows that the relative spectra have changed: 
the CIV broad absorption feature 16 years ago was wider in 
image D than in the other images, whereas it now looks the 
same, and the SilV+OIV] emission feature is even smaller in the 
VIMOS dataset (~l/3 compared to the average o f the other im- 
ages m easured from the continuum) than in the lAngonin et al] 
d 1 990 ) dataset (~l/2 compared to the average of the other im- 
age s). In addition, t here a ppears to be a slight slope difference in 
the lAngonin et al] ([1990) dataset when image D is compared to 
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the other images. However, this difference is not observed in the 
VIMOS dataset. Again some or all of these features could be ex- 
plained by three phenomena: galactic extinction, microlensing 
or intrinsic variability. We comment below upon the evidences 
against or in favor of such possibilities: 

- Galactic extinction : if there is a steeper slope in the con- 
tinuum of a quasar image when compared to the other im- 
ages, it can be inte rpreted as galactic extinction. In the 
Ango nin et al.l (Il990f) dataset this might be the case with 
image D, however, no slope difference can be seen seen in 
the VIMOS dataset. As galactic extinction is not expected 
to vary over time (especially on such a short time scale 
of 15 years) the slope should have persisted. Furthermore, 
iTurnshek et al.l (Il997l) showed that the reddenings between 
images B and C, as well as between images A and C were 
considerably greater than that between images D and C. 
Finally, differential extinction cannot explain the relative flux 
differences between the continuum and the emission lines 
when measured in different images. Therefore, we reject this 
interpretation. 

- Microlensing: the variation of the continuum slope of im- 
age D, together with the variation on the relative intensity 
of the emission lines in all four images (particularly image 
D) give a clear signal that microlensing is at work and is the 
best candidate for generating the observed changes. As ex- 
plained previously, microlensing can explain the chromatic 
continuum magnification. In this case, the image that suffers 
magnification of its continuum is also the one that showed 
a bluer spectrum, which is the most likely kind of chro- 
maticity induced by continuum microlensi ng. The differ- 
ences of the broad absorption features seen in lAngonin et al.l 
( 1990) were interpreted as microlensing from broad absorp- 
tion clouds. This suggestion has been tested with different 
microlensing simulations over t he years (e.g.. IHutsemekersI 
119931: iHutsemekers et al.lll994t ILewis & Bellelll998l) . sim- 
ulations in which the BAL profile of image D was repro- 
ducible. These differences are not seen any longer in the 
VIMOS dataset, which again is a result favoring microlens- 
ing in image D, as 15 years is roughly the expe cted timescale 
for m icrolensi ng variations in H 14 13+117 (IHutsemekersI 
1993). Finally, Oste nsen et al.l(ll997l) have shown, from their 
photometric monitoring between 1987 and 1995, almost par- 
allel light curves for the quasar images, except for image D. 
Image D exhibits a slightly higher amplitude in its overall 
variation, which, added to its spectral differences, makes a 
strong case for additional microlensing affecting image D. 

- Intrinsic variability : intrinsic variations of the quasar flux 
coupled with time-delays between its multiple images would 
also induce spectral differences between the images. Time- 
delays for this quasar have not been measured yet. Even 
though its expected value depends on the model, the sym- 
metric disposition of the images suggest s that it should be 
of the order of a month (e.g., model by IChae & Turnshekl 
1999). If the timescale of the quasar intrinsic variations is 
shorter than the time-delay for image D, then the brightness 
differences could also be explained by this phenomenon and 
should be seen in all four images at epochs separated by the 
time-delays. 

Although it is not possible today to rule out the quasar in- 
trinsic variability as the cause of the spectral differences be- 
tween the images (without doing a comprehensive and quanti- 
tative study based upon regular and precise spectrophotometric 



monitoring), we find it very likely that the spectral differences 
seen inH 1413+117 result from microlensing. Microlensing can 
explain all the feature differences seen in the spectra and is a 
natural explanation. 

7. Summary and conclusions 

We used the VIMOS -IFU device at the VLT to search for ob- 
jects in the vicinity of four lensed quasars for which a galaxy 
group or galaxy cluster was either confirmed or expected as part 
of the lensing system. Additionally, at least two of these lensed 
quasars were suspected to be under flux (de)magnification due to 
microlensing. We have described our data reduction techniques 
for integral field spectroscopy as an aid for future VIMOS users. 

We were able to measure flux ratio differences between the 
different quasar images in three of the systems: HE 0230-2130, 
RX J091 1.4+0551 and H 1413+117. Using the high quality 
spectra collected, our interpretation of the origin of these phe- 
nomena is microlensing in the three cases. Based on differences 
in the flux ratios between the different images when measured in 
the continuum or the emission lines, we conclude that: HE 0230- 
2130 is affected by microlensing magnification on the continuum 
of image D while image A (A1+A2+A3) of RX J091 1.4+0551 
is affected by a microlensing de-magnification of the continuum. 
Finally image D ofH 1413+117 shows a consistent history of 
having been affected by microlensing over the past 15 years. 
However, even though less likely as described in Sect. 6.3, these 
flux ratio differences could also be attributed to intrinsic vari- 
ability from the quasars coupled with time-delays between the 
images. 

In the three cases, the flux ratios between the quasar images 
show or have shown a non-zero slope, which is also consis- 
tent with microlensing. Nevertheless, in the RX J091 1.4+0551 
and HE 0230-2130 systems this phenomenon could also be 
explained by dust extinction due to the lensing galaxy or, 
even more likely, a combination of both phenomena. For 
H 1413+1 17, if there has been a slope difference between the im- 
ages, galactic extinction is not a likely explanation as this slope 
is not seen in the present VIMOS dataset. 

If, as we suspect, the three systems are affected by microlens- 
ing, future photometric and/or spectroscopic monitoring should 
fully characterize the phenomenon and allow derivation of the 
physical properties of the background sources. 

We have also detected absorption features from the two main 
lensing galaxies in the system HE 0230-2130 at z\ =0.521 and 
Z2=0.523. A third galaxy at similar redshift (z=0.518) has also 
been detected ~17" away from the quasar images. In the field 
around RX J09 11.4+0551 we have newly identified a galaxy, 
at a distance AO ~12" from the quasar images, and at redshift 
z=0.484, which makes it independent of the known galaxy group 
in the line-of-sight to this quasar. Even though not directly bound 
to the lensing galaxies, its vicinity to the quasar line-of-sight 
suggests that it might influence the total lensing potential. For 
H 1413+117 the higher resolution spectra clearly reveal the al- 
ready known absorption systems at redshifts ~1.44 and ~ 1 .66 
with high accuracy. 

Unfortunately, due to a combination of low exposure time, 
intrinsic PSF distortions and non-optimal observing conditions, 
the quality of the data obtained for B 1359+154 does not reveal 
any new information on the system or its surroundings. 
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